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ABSTRACT 

The nature of the Local Interstellar Cloud (LIC) is highly constrained by the com- 
bination of in situ heliospheric and line-of-sight data towards nearby stars. We present 
a new interpretation of the LIC components of the absorption line data towards e CMa, 
based on recent atomic data that include new rates for the Mg^ to Mg*^ dielectronic 
recombination rate, and using in situ measurements of the temperature and density of 
neutral helium inside of the heliosphere. With these data we are able to place interest- 
ing limits on the gas phase abundance of carbon in the LIC. If the C/S abundance ratio 
is solar, ~ 20, then no simultaneous solution exists for the A^(Mg I), A^(Mg II), A^(C II) 
and A^(C II*) data. The combined column density and in situ data favor an abundance 
ratio Ac /As = ^7^11 

We find that the most probable gas phase C abundance is in the range 400 to 800 
ppm with a lower limit of ~ 330. We speculate that such a supersolar abundance could 
have come to be present in the LIC via destruction of decoupled dust grains. Similar 
enhanced C/H ratios are seen in very low column density material, A^(H)< 10^^ cm~^, 
towards several nearby stars. 

Subject headings: ISM: clouds — ISM: abundances — ultraviolet: ISM 



1. Introduction 

The interstellar cloud that surrounds our solar system, the Local Interstellar Cloud (LIC), 
is a tiny, radius < 1 pc, low density interstellar cloud that has been detected in absorption lines 
towards many towards nearby stars (e.g., McClintock et al. 1978; Bruhweiler 1984; Genova et al. 
1990; Lallement & Ferlet 1997; Hebrard et al. 1999; Cry & Jenkins 2001; Wood & Linsky 1997). The 
cloud is partially ionized, ~ 20 — 40%, apparently due to the low intensity, relatively hard ionizing 
radiation field (Slavin &: Frisch 2002; Frisch & Slavin 2003, hereafter SF02). Direct observations of 
interstellar He^ which penetrates the heliosphere gives us tight constraints on its temperature and 
density: T = 6300 ± 340 K, n(He°) = 0.015 ± 0.0015 cm^^ (Witte 2004; Mobius et al. 2004). 

The star e CMa (HD 52089) is located in the downwind direction of the flow of ISM past 
the Sun {i = 239.8°, b = —11.3°, d = 132 pc). The absorption line data set for the line of sight 



^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS 83, Cambridge, MA 02138 
^University of Chicago, Department of Astronomy and Astrophysics, 5460 S. Elhs Avenue, Chicago, IL 60637 



- 2 - 



towards e CMa is particularly complete (Gry & Jenkins 2001, hereafter GJ). Using this dataset, in 
particular the ratios A'"(Mg II)/iV(Mg I) and N{C ir)/A^(C II), GJ put constraints on the electron 
density, n(e), and temperature of the LIC and of the other clouds seen in this line of sight. As we 
discuss, their method involved assuming a gas phase abundance ratio for C/S of 20, appropriate 
to Solar abundances. However, we show that when updated atomic rates are used for the Mg+ 
recombination coefficient, in particular for dielectronic recombination, this solution becomes no 
longer viable. This leads us to the conclusion that C has a siipcrsolar abundance in the LIC, 
Ac = 400 — 800 ppm (parts per million). This result, which is consistent with detailed models of 
the photoionization of the LIC (Slavin & Frisch 2006, hereafter SF06), then allows for a solution 
for n(e) which is consistent with all other data. High carbon abundances are also seen towards 
several other stars that sample the LIC. 



2. Data 

GJ presented data obtained with the GHRS on HST, and IMAPS on ORFEUS, of the ab- 
sorption lines from a number of ions observed towards e CMa, including Mg I, Mg II, C II, C II* 
and S II. Multiple velocity components are seen in this line of sight, but here we concentrate on 
the LIC component, referred to as component 1 by GJ, at a heliocentric velocity of ~ 17 km s^^. 
The derived column densities, along with the errors estimated by GJ, arc listed in Table 1. GJ 
determined a lower limit on A^(C II) directly from the observed absorption line, but needed to rely 
on an assumption for the abundance ratio of C to S to derive the upper limit on the column density 
because of the degree of saturation of the 1335 A C II line. Prom these data, using pre-1982 atomic 
rate constants, GJ obtained T = 5700 - 8200 K and n(e) = 0.08 - 0.17 cm-^. 

For the LIC we have the additional constraint given by in situ observations of the temperature 
of inflowing He'' determined from direct measurements inside of the solar system by the GAS 
instrument on the Ulysses space craft, T = 6300 ± 340 K Witte (2004, see also Mobius et al. 
2004). Because of the small column density of the LIC along the e CMa line of sight, A^(H I) ~ 
(3 — 5) X 10^^ cm~^, the gas temperature is not expected to vary greatly along the line of sight. 
In our photoionization models in which we use a realistic ionizing radiation field (SF02, SF06), 
the typical total variation in temperature is ~ 10%. We note that if GJ had required consistency 
with the LIC temperature in their analysis, their minimum electron density would have increased 
to n(e) ~ 0.11 cm~^ (up from 0.08 cm~^). 

3. Mg Ionization 

In ionization equilibrium the density ratio of Mg"'"/Mg° is determined by 

[r(Mg°) + C°'^(Mg°)n(H+)] n(Mg°) = a(Mg+) n(e) n(Mg+), (1) 
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where r(Mg'^) is the photoionization rate, C^'^iMg^) is the charge transfer ionization rate for 
Mg'' with H+, a(Mg^) is the total Mg"*" recombination coefficient (including both radiative and 
dielectronic recombination), and n(e) is the electron density. The dominant contribution to the 
photoionization rate for Mg'^ is the FUV background, which comes mostly from A and B stars. The 
spectrum and flux of the background has been estimated by a number of authors. We choose to 
use the results of Gondhalekar et al. (1980) because they make use of the TD-1 observations, which 
remain the best direct measurements of the background from stars and includes the whole sky. 
This is important because the sky is extremely anisotropic in the FUV and other observations have 
covered only brighter than average portions of the sky. Integrating the Gondhalekar et al. (1980) 
spectrum over the Mg'' photoionization cross section leads to a value of r(Mg°) = 4.5 x 10~^^ s~^. 
The uncertainty of the FUV spectrum is estimated by Gondhalekar et al. (1980) to be 25%. The 
Mg'^ photoionization cross section we use comes from Verner et al. (1996) and is based on results 
from the Opacity Project. Our value for the photoionization rate is significantly below the value 
6.1 X 10^^^ s^-*^ used by GJ, which was based on estimates by Jura (1974) and Mathis et al. (1983). 
The rate of charge transfer ionization of Mg*^ is taken from the fit provided in Kingdon & Fcrland 
(1996), which was based on data from Allan et al. (1988), and gives C*^'^(Mg°) = 6.5 x 10"^^ 
cm^ s-^ at T = 6300 K. 

The recombination rates, both dielectronic and radiative, are taken from new calculations. 
The radiative rates we use are recently calculated by Badnell^ using the AUTOSTRUCTURE 
code. These rates are considerably larger than those calculated by Aldrovandi &; Pequignot (1973). 
We use the dielectronic recombination rates calculated by this group as well (Altun et al. 2006), 
which are also larger than the older combined low temperature (Nussbaumer &: Storey 1986) and 
high temperature (Burgess 1965) rates. The resulting total Mg^ recombination rate for the LIC 
temperature range is Q;(Mg''') = 1.71[!^q 34 x 10~^^ cm^ s~^. As a result of the higher rates, the 
electron density inferred for a given observed ratio N(Mg II)/A'^(Mg I) is considerably lower than 
if the older rates had been used. It should be emphasized that there are considerable uncertainties 
in these theoretical dielectronic recombination rates. The low temperature (< 10^ K) dielectronic 
recombination rates in particular need to be confirmed by laboratory measurements and further 
theoretical calculations. Since we do not have any way of estimating the magnitude of the uncer- 
tainties in the rates we do not include the uncertainties in our numerical results and plots. Some 
comfort may be taken from the fact that other recently calculated rates (Gu 2004) are within ^ 30% 
of the Altun et al. (2006) results for temperatures of interest, though this does not guarantee their 
accuracy to that level. 

To use equation (1) to solve for the electron density, one needs to know n(H+). Based on our 

ionization models (SF02) and the observation that iV(H l)/iV(He I) > 12 (Dupuis et al. 1995) in 
the local ISM, we assume that n(II"^) = 0.9n(e). This should be good to about 10% and in any 
case the charge transfer term is small compared with the recombination term in the equation for 
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the electron density. 



4. C+ Fine Structure Excitation 

The ratio of C"*" ions in the excited J = 3/2 fine structure level of the ground state to that in 
the J = 1/2 level is another measure of n(e). The equilibrium populations of the two states are 
given by 

n{C+*) ^ 7i2n(e)+7i^n(H°) 

n(C+) A21 + 72in(e) + jgn{}P) ' ^ ^ 

where A21 = 2.29 x 10~^ is the radiative decay probability, and 712 and 721 are the excitation 
and de-excitation coefficients respectively for collisions with electrons and 7/3 and 721 are the 
corresponding coefficients for excitation by neutral H. The electron collisional excitation coefficient 
is 712 = ^'^^^1% ^' ^12 Gxp (— -fj^), where gi is the statistical weight (= 2) of the lower level, ^li2 is 
the collision strength and E12 is the excitation energy (see, e.g. Spitzer 1978). For partially ionized 
warm gas such as in the LIC the excitation of C+ by neutral H is small compared with excitation 
by electrons but for completeness we include it using rates from Keenan et al. (1986) and assuming 
a H'^ density of 0.2 cm~^, which should be roughly correct. We use the collision strength from 
Blum & Pradhan (1992), which is substantially smaller (^12 ~ 2.0 at T = 6300 K) than the Hayes 
& Nussbaumer (1984) value used by GJ (2.81). This smaller value is consistent within 5% with 
recent calculations by Wilson & Bell (2002). A smaller value for Q12 leads to a larger value for n(e) 
estimated from N{C ir)/N{C II). For T = 6300 K, 712 = 1.07 x lO"'^ cm^ 8"^ 



5. Results 

Following GJ, we determine n{e) from the overlap region that is consistent with both A^(C II*)/A^(C II) 
and iV(Mg I)/iV(Mg II). The updated atomic data we use leads to significantly different values 
from those found by GJ. Figure 1 illustrates the allowed region for n(e) as a function of T, based 
on equations (2) and (1) and the LIC temperature constraints. The allowed range of n{e) is much 
smaller than found previously primarily because of the limits on the temperature provided by the in 
situ observations. As can be seen from equation (2), the lower limit on n(e) from A'^(C II*)/A'^(C II) 
comes from the lower limit on the observed ratio. In the figure the dashed line, labeled "C II* /C II 
(C/S = 20)", uses an upper limit on N{C II) derived from taking the upper limit on N{S II) and 
an assumed C/S abundance ratio of 20. Clearly this value does not allow for a solution consistent 
with the Mg ion column densities and LIC temperature. The lower curve corresponds to a A^(C II) 
upper limit derived from assuming a C/S abundance ratio of 40, and allows a solution with the 
resulting limits of n(e) = 0.050 — 0.104 cm~^ when T is limited by the in situ He data. While C/S 
need not be as high as 40 to find an allowed region in the plot, a value of > 21 is required. 

An alternative method of placing constraints on the C/S abundance ratio is to use iV(Mg I), 
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N{Mg II) and T to derive n(e), and then N{C II*) and N{S II) to find ^cMs (^C = C/H and 
= S/H). Combining equations (1) and (2), we find 

Ac ^ ^2i+72in(e)iV(C IF) 

7i2n(e) N{S II) ^ ^ 

where n(e) is derived via equation (1). Observations with reasonable uncertainties exist for all 
of the quantities needed to determine Aq/As- A straightforward propagation of errors approach 
yields Aq/As = 46.8 ± 25.7 (including uncertainties in the temperature and its effect on the rate 
coefficients). The uncertainties involved in the expressions contribute in nonlinear ways, however, 
so we have in addition carried out a Monte Carlo simulation to estimate the errors. To do this we 
have used the 1-a reported errors for the observed quantities, assumed normally distributed errors 
and calculated the resulting Aq/As from equation 3 for many trials. The resulting probability 
distribution is very asymmetric with a long tail to high values of Ac/Ag. The result is Aq/Aq 
= 46.8^11'^, ■^here the error limits are those that contain 68.3% of the probability surrounding the 
peak in the distribution. The peak of the probability distribution function is at 37.2 and the mean 
(expected value) is 58.5. Alternatively we may state that the Ac /As > 21 with 95.45% confidence 
or Ac /As > 39 with 68.3% confidence (note that the long tail to high values causes this one- 
sided limit to exceed the value at the peak of the probability distribution). The high- value tail is a 
standard result for the ratio of two values in which the errors for each quantity are a relatively large 
fraction of the quantity. The inclusion of the uncertainty in the Mg^ photoionization rate broadens 
the uncertainty range in n{e) beyond that indicated in Figure 1. The Mg ionization calculations 
also yield the probability distribution for electron density for which we find n(e) = 0.064lQ Q3g, 
where again we report the limits that enclose 68.3% of the probability. For n(e) the peak is at 
0.052 cm~^ and the expected value is 0.067 cm~^. 

Determinations of the electron density for the LIC using independent sets of data, combined 
with models, lead Izmodenov et al. (2003), to find n(e) = 0.07 it 0.02. Similar calculations by 
Gloeckler &; Geiss (2004), but with somewhat different assumptions gave n{e) = 0.049 ± 0.016. 

Another approach for evaluating the C abundance, that is separate from the S abundance, 
uses only the more reliable A''(C II*) data. GJ placed limits on A''(H I) and iV(Htot) using several 
assumptions. Using the fact that the ionization of O and H are tightly coupled by charge transfer 
(Field &: Steigman 1971), observations of O I and an assumed gas phase abundance of 316 ppm, 
they find A^(H I) = 4.4^'^Qg x 10^^ cm~^. A recent determination of 0/H in the Local Bubble 
(Oliveira et al. 2005, also see Andre et al. (2003)) gives a value of 345 ± 19 ppm, which leads to a 
somewhat smaller total H I column density, A/"(H I) = 4.01q g x 10^^ cm~^. 

The ionization fraction of H in the LIC can be estimated in various ways. Using just the 
directly measured EUVE flux from stars, a lower limit of X(H"'") > 0.15 can be determined. This 
flux is inadequate to explain the high He ionization fraction relative to H found from nearby white 
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dwarf stars, however.^ Our models for the ionization of the LIC, which include the diffuse soft X-ray 

background as an ionization source and a contribution from radiation generated in an evaporative 
boundary between the LIC and the hot gas of the Local Bubble, predict ionization fractions of 
22 — 36%. For our best new models we find total H columns ranging from 4.1 — 5.5 x 10^''. While 
we cannot place firm upper limits on the H ion fraction, it appears unlikely that it significantly 
exceeds 40% unless some additional unknown ionizing source is present or the cloud is out of 
photoionization equilibrium (sec discussion below). Combining our H ionization fraction limits 
with the H I column determination leads to A^(Htot) ~ S.G^'^'g x 10^'' cm^^. This value yields a 
value for the S abundance of Ag 15 ppm, with 50% uncertainties, in agreement with solar values 
of 12.3 - 15.5 (Asplund ct al. 2005, AGS). 

We can express the C abundance as a function of the observed value of N{C 11*) as, 

iV(Cir) A21 f Ar(Cir) \ /5 X 10i^cm-2\ /o.08cm-3\ 



iV(Htot) 7i2n(e) Vl-3 x 10^2 cm'V V A^(Htot) J\ n{e) 

where we have left out the less important terms for collisional de-excitation and excitation by H°. 
This expression makes it clear that the C abundance in the LIC is quite high unless N{C, II*) is at 
the low end of its range while A^(Htot) and n(e) arc the high ends of theirs. In this expression n(e) 
should be the line of sight averaged value, which is highest at the cloud edge and decreases inward 
on the line of sight. For our models we find that the expression leads to about an 8% overestimate 
for Ac when n(e) is interpreted as the electron density at the position of the heliosphere. If we 
take an upper limit on n(e) = 0.09 cm~^, consistent with our derivation above, and an upper limit 
of iV(Htot) = 8 X 10^^ cm-2, and the lower limit of iV(C II*) = 1.1 x 10^^ ^m-^ we get 334 ppm 
as a lower limit on the C abundance. This lower limit is consistent with the solar abundance from 
Grcvcssc & Sauval (1998), but not with the more recent results of AGS. Using a value of A^(Htot) 
as high as 8 X 10^ implies a S abundance of ^ 13 ppm to be consistent with the observed upper 
limit for A''(S II). Sulfur is expected to have little to no depletion in warm clouds (see e.g., Welty 
et al. 1999), and a solar abundance of 13.8 ppm is found by AGS. 



6. Discussion 

An implicit assumption in our analysis for the electron density is photoionization equilibrium. 

The timescale for Mg I ionization is very short, < 10^ years, while the recombination timescale 
for Mg II is considerably longer 3 x 10^ yr. The most likely scenario for the cloud to be out of 
photoionization equilibrium is that it is recombining and cooling from a previously heated (e.g., 
shocked) state (see Lyu & Bruhweiler 1996). In that case the electron density we infer would be an 
underestimate. This would both lower the inferred value of Ac/Aq and the C abundance. However, 
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0.38 respectively (Wolff et al. 1999, and references therein). 
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Jenkins et al. (2000) argue, based on the relative column densities of Ar I and H I in the LIC, that 

photoionization dominates the ionization rather than the collisional processes that should dominate 
in cooling and rccombining gas. It remains to be seen whether or not a non-equilibrium model could 
be found that would be consistent with the Ar I data as well as all the other column density data 
for the LIC and the in situ data on He°. We intend to explore this possibility in the future. 

Observations of low column density, A'^(H)< 10^^ cm~^, ISM towards other nearby stars, are 
consistent with our conclusion of supersolar C abundances. Using N I+N II as a proxy for the total 
H column and the AGS Solar abundance for N for the lines of sight towards Capella, WD 1254+223, 
WD 1314+293, and WD 1634-573 yields Ac ~ 340 ppm, with a range 83 - 550 (see data in Wood 
et al. 2002; Lehner et al. 2003; Oliveira et al. 2005; Kruk et al. 2002). Three of these stars show 
Ac > 333 ppm. 

7. Conclusions 

We have shown that using updated atomic data and the in situ data on the temperature of 
the gas in the LIC, we can put interesting constraints on the C abundance in the LIC. Unless 
a number of observed or inferred data values are at the far ends of their uncertainty ranges, it 
appears that C is not only undepleted in the LIC, but has an abundance of Ac ~ 400 — 800 ppm, 
considerably larger than the solar value. In addition, the ratio of C to S abundance appears to 
be necessarily larger than the solar ratio of ~ 20 with values in the range of 40 — 50 favored by 
the constraints on n(e) in the LIC. We speculate that this supersolar abundance of C, if it truly 
exists, may be indicative of a local enhancement in dust in the past. In this picture the gas phase C 
abundance was enhanced by the passage of an interstellar shock that destroyed the carbon bearing 
grains more effectively than other, e.g. silicate, grains. This could be related to the unusual dust 
size distribution in the LIC as indicated by Ulysses observations of large grains penetrating into 
the heliosphere (Gruen et al. 1994; Prisch et al. 1999). 

Further investigations into the C abundance in the LIC should be undertaken to confirm 
our findings, especially studies of other lines of sight that show the LIC velocity component in 
the absorption lines of C, S, Mg and other elements. The LIC is the best observed low density 
interstellar cloud and thus can provide us with valuable information on this phase of the interstellar 
medium. The C abundance and its evolution in the warm diffuse ISM has many implications for 
the nature of interstellar dust and metallicity in the Galaxy. 

This research was supported by NASA grants NAG5-12571 and NAG5-13107. The authors 
thank Daniel Savin and Gary Ferland for very helpful discussions and the anonymous referee who 
helped improve the clarity of the presentation. 
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Fig. 1. — Constraints on the electron density and temperature in the LIC based on column density 
measurements towards e CMa and in situ observations of the temperature. The thick lines show the 
constraints from the observations of Mg I and Mg II. The upper, thin solid line is the upper limit 
on the electron density using the lower limit on iV(C II) derived directly from the observations. 
The dashed and dashed-dot lines show the lower limit on n{e) that result from an upper limit on 
N{C II) under the assumptions of a C/S abundance ratio of 20 and 40 respectively. The vertical 
dotted lines indicate the temperature limits based on in situ observations of He°. 



- 11 - 



Table 1. Observational Constraints 



Observed 


Observed 


Reference 


Quantity 


Value 




iV(C II) (cm-2) 


1.4 - 2.1 X 10^'^ 


1 


N{C II*) (cm-2) 


1.3 ± 0.2 X 10^2 


1 


N{0 I) (cm-2) 


i.4i°:^ X 10^^ 


1 


A^(Mgl) (cm-2) 


7±2 X 10^ 


1 


A^(Mg II) (cm-2) 


3.1 ±0.1 X 10^2 


1 


A^(S II) (cm-2) 


8.6 ±2.1 X 10^2 


1 


r(K) 


6300 ± 340 


2 


n(He I) (cm-3) 


0.015 ± 0.0015 


2 



References. — (1) Gry & Jenkins (2001) (Values 
shown are for the LIC component towards e CMa, 
component 1.), (2) Witte (2004). 



